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SUMMARY 

As a part of the dynamical heating of the upper at- 
mosphere, which has been emphasized by Hines, the 
acoustic heating of the upper air is reinvestigated, par- 
ticularly the warming of the polar night mesosphere. 

It is shown that because of the relatively small trans- 
missivity of the thermosphere for long period sound waves, 
most acoustic heating does not take place above the 200 km 
level, but rather below 100 km or around the mesopause, 
and that although the acoustic output power from the polar 
jet stream is more than 10 times larger in winter than in 
summer, acoustic heating is not sufficient to compensate 
for the cooling rate around the mesopause in winter, unless 
the wind velocity of the polar night jet stream continually 
exceeds 200 m/sec. 
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ACOUSTIC HEATING OF 
THE POLAR NIGHT MESOSPHERE 


by 

Kaichi Maeda 

Goddard Space Flight Center 


INTRODUCTION 

An effective heating mechanism for the polar mesosphere, which is warmer in winter than in sum- 
mer, is Kellogg 1 s chemical heating through the recombination of atomic oxygen (References 1, 2, and 
3). Although the velocity of subsidence required for this process is significantly less than that re- 
quired for adiabatic heating alone, it must still be of the order of 0.1 cm/sec near the 60-90 km level. 

On the other hand, theoretical calculations by Haurwitz (Reference 4) indicate an upward rather 
than a downward flow in this region. And although Sawada (1961) has questioned Haurwitz’s omission 
of certain inertia terms in his basic equations, a recent analysis by Kochanski (Reference 5) of 
rocket data relating to mesospheric circulation shows good agreement with Haurwitz's results. 

Rocket soundings by Nordberg and Smith (Reference 6) have occasionally measured very warm 
wintertime mesosphere temperatures even in the middle latitudes, where subsidences are unlikely 
but very strong westerlies in the 40-80 km region appear to be associated. 

Hines (References 7 and 8) has suggested that the increase of temperature in the polar night 
mesosphere could be explained by dynamical heating. Gassard (Reference 9) and Blamont and Jage 
(Reference 10) have found evidence of upward propagation of tropospheric disturbances into the 
ionosphere. 

The possibility of acoustic heating of the ionosphere was discussed by Daniel (Reference 11), who 
considered the noise produced by the sea surface as an energy source. However, Eckart (Reference 
12) showed that pressure waves from this source are completely ineffective in transferring energy up 
to such high altitudes. This has recently been confirmed by Cook and Young (Reference 13). The 
amount of acoustic heating of the ionosphere is also estimated by Golitsyn (Reference 14), using 
tropospheric sources. He arrives at an upward energy flux of the order of 10" 2 erg/ cm 2 sec, which 
is insufficient to cause the observed heating. Moreover, Golitsyn did not consider the reflection of 
acoustic energy due to the vertical temperature distribution in the atmosphere. 

In this paper, the contribution of aerodynamically produced atmospheric acoustic waves to the 
warming of the polar night mesosphere is reexamined, with the polar night jet stream in the 
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stratosphere as the source, and with the thermal structure of the polar upper atmosphere, playing a 
part in propagation and absorption. 

It should be noticed that the largest portion of the kinetic energy of the atmosphere is in the 
troposphere, where the air density is greatest, and is of the order of 10 3 ergs/cm 3 . If this tropo- 
spheric kinetic energy is not trapped effectively within the troposphere, the earth’s atmosphere at 
ionospheric heights should reach very high temperatures of the order of 10 5 °K. However, Charney 
and Drazin (Reference 15) have shown that the escape of large amounts of planetary wave energy from 
the tropospheric layer into the upper atmosphere is prevented throughout most of the year by the large 
westerly zonal wind systems above the tropopause, except for short periods in the spring and autumn. 


KINETIC ENERGY DENSITY OF THE POLAR STRATOSPHERE 

Owing to the development of the polar night vortex and jet stream, the kinetic energy of the polar 
stratosphere increases remarkably during the polar nights. As may be the case for the solar corona, 
some part of this kinetic energy may be transferred upward by means of internal pressure waves. 



1957 1958 

DATE ( months ) 

Figure 1— Monthly averages of kinetic energy density 
above Ft. Churchill (59°N geographic latitude) at 30 
mb (solid lines) and at 50 mb (dashed lines) during the 
period from July 1957 to July 1958. 


According to Boville, Wilson and Hare 
(Reference 16), the kinetic energy density of zonal 
wind systems between 60°N and 80^ at 25 and 500 
mb increases significantly during the wintertime 
(from December until March). The kinetic energy 
density corresponding to the winter peak of the 
basic zonal flow is of the order of 200 and 10 3 
ergs/cm 3 at the 25 and 500 mb levels, respectively. 

Figure 1 shows monthly averages (Refer- 
ence 17) of kinetic energy densities during the 
period from July 1957 to July 1958 at 30 mb 
(solid lines) and at 50 mb (dashed lines) above 
Ft. Churchill (59 °N geographic). This increase 
of kinetic energy in the upper stratosphere in 
wintertime is mainly due to the development of 
the polar night jet stream (Reference 18). 

The dissipation of the energy of jet streams, 
as well as any strong wind system in the at- 
mosphere, is mainly due to friction in the bound- 
ary layers, which cause the turbulent motion to 
decay into small eddies in the surrounding 
mediums. These eddies decay further into ther- 
mal motion and cause local heating. However, 
some part of the kinetic energy in turbulent mo- 
tion in jet flow can propagate into the upper 
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atmosphere in the form of pressure waves. Owing to the exponential decrease of air density with 
height, the amplitudes of those propagating pressure waves grow rapidly at high altitudes and ap- 
proach shock- wave magnitudes, i.e. dissipation becomes non-linear. 


THE EQUATION OF INTERNAL ATMOSPHERIC WAVES 
Fundamental Equations 

The equation of internal atmospheric waves is derived from the following three fundamental 
equations: the equation of motion 


Pv 

Dt 



+ g 


f ; 


the equation of continuity 


jo 

Dt 




and the equation of thermodynamics 


Dt 



(i) 


( 2 ) 


( 3 ) 


Since the source of disturbance is extended along a certain latitude, we can treat the problem in 
two dimensions and specify the coordinates within a meridional plane. Let x be the horizontal (south- 
ward) and z the vertical (upward). 

In the first approximation, the Coriolis force and external forces, except gravity, are neglected. 
The equation of motion (Equation 1) then becomes 

du _ _ dj) 

P 0 d t <3 X 


and 


The equation of continuity (Equation 


dw _ djg 
Po dt d z * 


2) for the first order approximation is 


dp dp o 

dt + w TT - - Po* ’ 


( 4 ) 
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where X V(x, z,t) is the velocity divergence as defined in the previous section. Using Equation 4, 
the equation of thermodynamics (Equation 3) can be written as 

dj) f -y \ 

dt + Ic'x-wgJ p Q - 0 

The waves travelling in the flat atmosphere can be found by assuming that u , w , p and p are pro- 
portional to e i <“ t “ kx >. Then, the equation for vertical change of x(^,z) is given by eliminating u , 
w, p and p: 


d 2 X _1 

dz 2 c 2 

Since c 2 = >gH, this can be written as 


/dc 2 \ 


'to* 

,9 g 2 k 2 / dc 2 \ 

[~57 -rgj 

dz + 

_c 2 

k T a) 2 c 2 1 dz + (T-l)eJj 


d 2 y 1 dv 

-772 + ¥(»'- i ) af + 


.^' k2+ l^¥ ( H ' + x -r). 


o , 


where H' - dE'dz. 


The Diagnostic Diagram 

The wave equation (Equation 5) can be written 


(5) 


dz 


f - 2N ^7 + M 2 x = 0 


( 6 ) 


where 


N 


1 




) . 


M 2 




The solution of the differential equation (Equation 6) is given approximately by 


z) - e N * ( Ae nz + Be +nz ) , 

where A and B are constant and can be determined by boundary conditions. 


( 7 ) 

( 8 ) 

(9) 
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It should be noted that the solution (Equation 9) of Equations 5 or 6 can be classified into cellular 
and non-cellular types*, depending upon the relative magnitudes of n andM, i.e., 

(i) cellular solution 


n - i?j , t ) 2 


M 2 - N 2 for M 2 > N 2 , 


(10) 


(ii) non-cellular solution 


n 2 = N 2 - M 2 for N 2 < M 2 . (H) 

The cellular waves have propagation vectors in both the horizontal (k) and vertical (^) direc- 
tions, but non-cellular wave has no propagation vector in the vertical direction and simply propagates 
in the horizontal direction. 

Since M is a function of the horizontal wave number, k, Equation 10 gives the relation between two 
wave numbers in the vertical (77) and the horizontal direction (k): Equation 10 is written as 

T j 2 C 2 k 2 c 2 ^ 2 -^ B 2 ) 

" 2 ""a 2 1 (12) 


where 


“a 2 - (* _ ji ar ) or "a - ( 13 ) 

“b 1 = ^ (tt + (y~ ^e) or u b = r + " 7) : 04) 

co A and are called the critical acoustic angular frequency and the Brunt- Vaisala angular frequency. 
These are both functions of the atmospheric temperature distribution in the vertical direction and are 
related to the stability of layers in the atmosphere. 

Since 6 j a > co B we can see from Equation 12 that, if <u > o> A , the propagation surface is an ellipsoid 
while if co < co B , it is an hyperboloid. In both cases, the rotation axis (z-axis) is vertical (Figure 2). 
The waves corresponding to co > co A are called sonic or acoustic waves and the ones for < co B are called 
thermobaric waves (Reference 19) or internal gravity waves (Reference 20) respectively. 

Since 77 2 becomes negative for co B < w < co k , no cellular wave exists for this frequency range. This 
situation can be seen from Figure 3, in which are plotted as a function of k the two curves of w, cor- 
responding to t } 2 - n 2 - 0 . Similar curves are shown in Figure 4, where the wave period (t in sec) is 

•The nomenclature of cellular and non-cellular is taken from Martyn’s paper (Reference 21) and this corresponds to internal and external 
waves discussed by Charney and Drazin (Reference 20). 
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Figure 2— Propagation surfaces (I in es) of interna! atmospheric waves in the isothermal 
atmosphere with scale height 8 km (T = 0°K), given by Equation 12; the solid and 
dashed lines correspond to sonic waves (acoustic waves) and thermobaric waves (in- 
ternal gravity waves), respectively. 


plotted against the horizontal wave length instead 
of co and k, respectively. r A and r B correspond 
to anduJ B , respectively. 

Figures 3 and 4 are called diagnostic dia- 
grams of atmospheric waves (Reference 21), in 
which a domain indicated by n 2 > 0 corresponds 
to the non- cellular mode while two domains 
n 2 < 0 correspond to the cellular modes. 

APPROXIMATE TRANSMISSIVITY OF 
THE POLAR ATMOSPHERE 

Since the non- cellular (or external) waves, 
the period of which is limited between r A and r B , 
have no vertical propagation vector, these waves 
neither contribute to the vertical energy 
transfer nor exist locally in free atmosphere. 



Figure 3— Diagnostic diagram of interna! atmospheric 
waves in the isothermal atmosphere with scale height 
8 km (T = 0°K). 
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Figure 4—DiagnosHc diagram of internal atmospheric waves in the isothermal at- 
mosphere with scale height 8 km (T = 0°K). 


Among two kinds of cellular (or internal) waves, on the other hand, the thermobaric (or internal 
gravity) mode can propagate with larger periods than the Brunt period r B , while the sonic (or acoustic) 
modes can exist only below a period r A . 

As can be seen from Equations 13 and 14, these two critical periods are a function of temperature 
and vertical temperature gradients in the atmosphere (i.e., the lapse rate of atmospheric tempera- 
ture). Therefore, these values differ not only with height but also with season. In Figures 5 and 6, 
r A and t b are plotted against altitude for polar atmospheres above Fort Churchill in summer and in 
winter, respectively. 

The atmospheric models of the two seasons are essentially the same as the ones reported by 
Stroud et al. (Reference 22; see also Figure 13, Reference 23 and Figure 1, Reference 2). In a 
domain between two critical periods, i.e., t a (solid curves) and r B (dashed curves), only the non* 
cellular mode can exist. In other words, this is a forbidden region for both acoustic and thermobaric 
modes of internal atmospheric waves. Therefore, we can see that no thermobaric waves with periods 
shorter than 370 sec can propagate across the 50-km level and that no sonic waves with periods 
longer than 250 sec can propagate through the layer around 80 km in summertime. But in winter both 
waves are able to propagate through these layers and the critical period for each approaches 320 sec 
for thermobaric waves and 290 sec for sonic waves, respectively. 


To see the transparency (or transmissivity) of the atmosphere for internal atmospheric waves 
more quantitatively, the wave equation (Equation 5 or 6) should be transformed into a more simple 


7 





PERIOD (sec) 


Figure 5— Critical acoustic period, r A , and Brunt- 
Vaisala period, r B , as a function of altitude for the 
summer polar atmosphere. 



Figure 6— Critical acoustic period, t a , and Brunt- 
Vdisala period, r B , as a function of altitude for the 
winter polar atmosphere. 


form. By substituting 


y 


£ 


dz 

H 


and 


X<*) = c^ 3 0(y) , 


the term d x /dz in the wave equation can be eliminated; then Equation 6 can be written as 


dV(y) 

dy 2 


- K 2 4>{y) 


(15) 


where 


K 2 = 7j 2 • H 2 . 


As can be seen from Equations 7, 8, and 10, rf is a function of cr, k, H, and H' ( - dH'dz). 
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The one-dimensional wave equation (Equation 15) indicates that, if K 2 > 0, the solution 4>( y) has a 
wave nature and the corresponding waves can propagate along the y-axis (z-axis vertical direction) 
within the domain, where this condition is satisfied. On the other hand, if K 2 £0 (Kis imaginary) i.e., 
no internal wave exists and vertical propagation of waves is inhibited by the barrier, where K 2 be- 
comes negative. Thus, the propagation of internal waves is limited within regions where K is real. 

For the acoustic mode of atmospheric waves with periods 100 sec and 10 sec, K 2 is plotted against 
z in Figures 7a and 7b for summertime and in Figures 8a and 8b for wintertime, respectively. In 
these figures, the horizontal wavelength is taken as a parameter. 

Since \ = 2n/ k, where k is a wave number in the horizontal direction, there is a shortest (critical) 
horizontal wavelength k c for a given period of the waves, corresponding to a largest horizontal wave 
number k c . From Figures 7 and 8, we can see that the waves propagating close to the vertical axis 
can propagate upwards as far as r < r A , while the waves propagating in an oblique direction meet a 
barrier where K 2 becomes negative even if the condition r < r A is satisfied. If there are more than 
two barriers, those obliquely propagating waves are trapped between these barriers, forming the 
"duct" in the horizontal direction. It should be noted from Figure 2 that k = 0 (i.e. k - <*■) represents 
the wave propagation in the vertical direction. 

Using the diagram shown in Figure 2 for the conditions at each level, the critical angle 6 C can 
be drawn as a function of altitude for a given period, r, of the waves, provided that r < r . The angle 
& c is defined such that the waves propagating obliquely with zenith angle larger than Q c cannot enter 
into the above layer, simply because K 2 becomes negative beyond this angle. In other words, 8 c can 
be regarded as the aperture of the allowed cone around the vertical axis at each level in the atmos- 
phere. Therefore, the relative transmissivity, f t (r, z), of each layer for vertical propagation of 
the atmospheric acoustic wave is approximately given by (# C /V 2 ) 2 > provided that is expressed 
in radians. In Figure 9, f t (r, z ) is shown for the wave with a 100 sec period normalized to the con- 
dition at 25 km. Similar curves are obtainable for any waves with r < r A . From Figure 9, we can 
see that the atmosphere above mesopause, which is nearly 80 km in polar regions, is always opaque 
for atmospheric acoustic waves. 

It should be noted that the wind system in the higher layers does not effect f t (r, z ) very much, 
because an oblique axis tilted windward direction can be regarded as the direction of maximum 
propagation. 

UPWARD ENERGY TRANSFER BY PRESSURE WAVES 
Acoustic Waves Generated Aerodynamically in the 


Sound waves generated aerodynamically are distinct from sound produced by the vibration of solids. 
A strong jet flow in the free atmosphere produces sound waves due to turbulent motions in the bound- 
ary layers between the core of the jet and the surrounding air. Lighthill (References 24 and 25) in- 
vestigated this problem and derived a theoretical expression of the power of sound waves as a function 
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of the jet stream velocity. His theory was 
further extended by Proudman (Reference 26) 
to the generation of sonic waves by isotropic 
turbulence. 

Although the acoustic heating of the solar 
corona is discussed by Schatzman (Reference 
27), similar calculations cannot be applied to 
the earth’s atmosphere, because the situation 
there is entirely different from that of solar 
atmosphere. In the earth’s upper atmosphere, 
the effect of magnetic fields on the generation 
of sound noise by isotropic turbulence (Refer- 
ence 28) is important as well as their effect on 
the dissipation processes as in the case of 
solar corona. 

Since the aerodynamical production of 
acoustic waves has been developed in detail by 
Lighthill and Proudman, only their results will 
be given, and, in turn, applied to the present 
problem in the earth’s upper atmosphere. 

The acoustic power output per unit mass 
of air P, in ergs/gm-sec, is given by 

P - aeM 5 for M « 1 (16) 



Figure 9 — Relative transmissivity of the polar atmosphere, 
f t (t, z) for atmospheric acoustic waves with a period 
t = 100 sec. The solid and dashed lines correspond to 
wintertime and summertime, respectively. 


where a is a nondimensional constant, m is Mach number e is the mean rate of energy dissipation per 
unit mass. 

According to Batchelor (Reference 29, p. 103), e is given by 


e 




where A is a nondimensional constant, which is of the order of unity— about 1.8, l a characteristic 
length which can be taken as the reciprocal of the wave number at which the maximum of the energy 
spectrum occurs, and u the root mean square velocity of the turbulence. 
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If the "double correlation function" given by Heisenberg (Reference 30) is correct, the constant 
a in Equation 16 has the value 38. Therefore Equation 16 is 

u 3 

P = 70 — * (M) 5 


where M = u/c, u and c are the mean velocity of jet flow and the sound velocity respectively. 

In the upper atmosphere, where the Reynolds number is very large, the spectrum of turbulence 
consists generally of the following ranges: 

1. The energy-producing range, which contains large eddies, 

2. The inertial subrange in which energy is neither created nor destroyed; and 

3. The energy-dissipating range, which contains small eddies. 

As pointed out by Proudman (Reference 26), the turbulence which makes an appreciable contribu- 
tion to the generation of sonic waves at large Reynolds numbers belongs to the eddies that do not dis- 
sipate energy. In other words, the main source of aerodynamically produced sonic waves is turbulent 
motion around the jet flow; this turbulence belongs to (1) energy producing range, and to (2) inertial 
sub-range. Therefore, following Batchelor (Reference 29, Chap. 7), we can assume 

io 3 s i & io 4 

and 

10 ^ u £ 3 

Distribution of turbulence along the jet stream can be assumed uniform, then the acoustic output 
per unit area along the jet (in erg/cm 2 -sec) will be 

7 • 10” 6 u 3 M 5 p Q < I(u) < 3 * 10" 3 u 3 M 5 p Q , 

where /? 0 is the mean density of air around the source. In Equation 17, u is in cm/sec. The acoustic 
output I (u ) is shown in Figure 10, assuming the height of the source to be 25 km. 

From this figure we can see that the production of sonic waves decreases drastically below a 
wind velocity of the order of 50 m/sec. Since the derivation of Equation 17 is based on the as- 
sumption that Mach number is small (Reference 25), the output power beyond u = 150 m/sec might 
be incorrect. 


13 



Attenuation of Atmospheric Waves 


10,000 


Any atmospheric wave is attenuated in air by the fol- 
lowing effects: 


1. friction due to the viscosity of air, 

2. thermal conduction of air, 

3. radiation through the atmosphere, 

4. dispersion due to the inhomogenity of air in the 
atmosphere. 

As was shown by Rayleigh (Reference 31), effects of ra- 
diation and dispersion are generally small for the acoustic 
waves as compared with the effects of friction and thermal 
conduction. 


Consideration of these effects leads to the so-called 
classical Stokes-Kirchhoff formula for the attenuation coef- 
ficient a (r) (Reference 31), which is given by 


a(r) 




(18) 


where a 2 and v are the coefficients of thermal conductivity 
and kinematic viscosity of air; 


1,000 b 



u ( m/sec ) 

Figure 10 — Acoustic output power of the 
polar stratosphere, !{u) , plotted against the 
mean wind velocity of the jet stream, u , 


a 


2 


C p 
P H 


and 


The heat conductivity of air, k , is the order of 6 x 10" 6 cal/cm-sec°K at NTP; and the specific 
heat C p is nearly 0.27 cal/gm'TC below the mesopause and increases with height to 0.3 cal/gm°K above 
350 km approximately (Reference 32). 

The viscosity of air, p, is of the order of 1.7 x 10“ 4 poise at NTP. Using these numerical values, 
the attenuation coefficient (Equation 18) becomes 


a(r) 


2.4 x 10" 16 


pr 


2 
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As was shown by Golitsyn (Refer- 
ence 33), the long period pressure 
wave undergoes an additional dissipa- 
tion owing to the Joule-loss of the 
current induced by fluctuations of the 
conducting medium of ionosphere in 
the earth's magnetic field. 

In Figure 11, the attenuation 
factor 


f a ( T » Z > 


exp - 



z * ) dz 


(19) 


is shown for four different periods of 
waves (1, 10, 30 and 100 sec). The 
dashed lines in Figure 11 are taken 
from Golitsyn’s results, which took 
account of the Joule loss (Refer- 
ence 33). 

An estimate of the rate of energy 
dissipation per unit time per unit 
thickness per unit cross section in 
the upper atmosphere, dE t/ /dz in 
erg/cm 3 -sec, is made by the follow- 
ing approximation: 



100 200 300 400 

ALTITUDE ( km ) 

Figure 1 1 —Attenuation factor, f (r, z) , of atmospheric acoustic 
wave versus altitude, z for r = 1, 10, 30 and 100 sec. The solid 
lines are calculated by the classical formula and the dashed lines 
are given by Golitsyn, with the Joule loss in the ionosphere taken 
into account. 


fit 

dz 



i(ai) f t (a;, z) dz f a (a), z) d^ , 


( 20 ) 


where 


f (oj, z), the attenuation factor, is given by Equation 19 and shown in Figure 11; 

f t (co 9 z) is the relative transmissivity of atmosphere for the wave with angular frequency w 
at the altitude z, which is shown, for example, in Figure 9 for r = (2n/o)) = 100 sec; and 

i(o;) is the differential frequency spectrum of the wave, which will be assumed to be 


i (uj) - i 0 oT m m > l* 


(21) 
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From Figures 5 and 6, the lower limit co A of the integral in Equation 20 can be taken as follows: 


(i) summer 


"A - "A ( z ) 

for 

0 1 

z <20 

km, 

0.023 

for 

20 1 

Z <70 

km, 

“a ( z ) 

for 

70 £ 

z <85 

km, 

0.025 

for 


z >85 

km; 

(ii) winter 





"A = "A ( z ) 

for 

0 < 

z <30 

km, 

0.024 

for 


z >30 

km. 


From this figure, we can see that the change of co A with height is small. Thus, co A is taken as a con- 
stant value (0.024) for both seasons. 

Strictly speaking, the right hand side of Equation 20 should be multiplied by a geometrical factor, 
which is inversely proportional to the distance from the source. The width of the source— the polar 
night jet stream— is more than several hundred kilometers (Reference 18) and the boundary is widely 
diffused horizontally. In this respect, the geometric attenuation factor is not important and can be 
neglected in the present approximation. 

Since the waves are propagating upwards, f t O, z) in the integrand of Equation 20 may be taken 
as a constant value above the minimum transmissivity level. However, the formation of wave-ducts is 
highly probable above this or between those levels of minimum transmissivity, for example, between 
30 and 70 km in summer (see Figure 5). As can be seen from the result, the choice of f t (a >, z) 
above the minimum transmissivity level (if it is constant or variable as shown in Figure 9, for ex- 
ample), does not affect the present calculation very much. 

A constant i 0 in Equation 21 is given by 


I 


o 



i Q of m dcr 


co 


1 

A 


m - 1 


( 22 ) 


where I 0 = l 0 (u) is the output power of the acoustic wave at the level of the source, z 0 , given by 
Equation 17 and shown in Figure 10. 

In Figure 12, the relative dissipation rate of energy (dE t /dz) /l 0 is plotted against altitude z for 
two cases of the spectrum, m = 2 (heavy lines) and m = 5 (light lines). The solid and dashed lines cor- 
respond to the wintertime and the summertime polar atmosphere, respectively. The two curves in 
the upper portion of Figure 12 are calculated by neglecting the transmissivity of the atmosphere, i.e. 
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Figure 12— Relative dissipation rate of energy flux transferred by atmospheric acoustic waves, where 
*0 = * o ' s the acoustic outputat the source as shown in Figure 10. Solid and dashed lines correspond 
to wintertime and summertime polar atmospheres, respectively; heavy and light lines stand for the 
hard spectrum ( m = 2) and soft spectrum ( m - 5) of the waves, respectively. 


f t ( go , z) is assumed to be unity. The latter corresponds to the estimation made by Golitsyn (Ref- 
erence 14). 

If all those dissipated energies of acoustic pressure waves are assumed to be converted into 
thermal energy of the air, the rate of heating at each level, dT(z)/dt in °C/sec, is given by 


dT 

dt 


1 

C P P( z ) 



d^ 

dz 


(23) 


where J 0 is Joule’s mechanical equivalent of heat, 4.2 x 10 7 ergs/cal; C p is the specific heat of air at 
constant pressure, -0.3 cal/gm°C; and p( z) is the mean atmospheric density at the altitude z. 

The results are shown in Figure 13, which corresponds to Figure 12 for a unit output power 
source. Since the latter is shown in Figure 10, the following points may be observed: 

If the wind velocity around 25 km level is of the order of 30 m/sec, as in the usual summertime 
polar stratosphere, we can see that the acoustic power at the source is less than the order of 10 5 
erg/cm 2 -sec (from Figure 10). Then rate of maximum heating, dT/dz is, from Figure 13, of the order 
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icr 10 icr 9 icr 8 io -7 io -6 icr 5 icr 4 io ~ 3 


J_dT / °C \ 

Io d z \ sec - I 0 J 


Figure 13— The rate of atmospheric 
heating in each layer due to the 
dissipation of acoustic energy flux 
*n (erg/cm 2 sec) versus altitude z. 
The lines are defined as in Fig- 
ure 12. 


of 10 8 x 10 5 °C/sec or 10 8 °C/day at the maximum, which height is around 75 km if the wave 
spectrum is hard ( m = 2), or 100 km if the spectrum is soft (m = 5), respectively. 

Oh the other hand, if the wind velocity around 25 km exceeds 100 m/sec, which is steadily ob- 
served in the polar night jet stream (References 17 and 34), the output power I 0 is somewhere be- 
tween 10 3 and 1 erg/cm 2 sec (Figure 10). The corresponding maximum warming rate is, then 
2 x (IO -5 to 10~ 8 ) °C/ sec or 2 x (1 to 10“ 3 ) °C/day at the altitude of 75 km for m = 2, or around 100 
km for m = 5 . 

If the transmissivity of the atmosphere is disregarded, as in Golitsyn’s estimation, then from the 
two upper curves shown in Figure 13, the rate of maximum heating is of the order of 10“ 10 °C/sec 
(10~ 5 °C/day) in summer, and 10 -2 to 10~ 5 (10 3 to l°C/day) in winter. The height of maximum heat- 
ing is around 260 to 300 km which corresponds to m = 2 to 5. 


CONCLUSION 

Although the estimated acoustic energy flux available for upper air heating has a wide range as 
shown in Figure 10 — due to ambiguities of several parameters such as the mean velocity of turbulent 
motion around the jet stream, u; the effective size of eddy in this turbulence, l; and the frequency 
spectrum of the turbulence, the following remarks can be made from the present calculation. 
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Because of the temperature gradient, the atmosphere above the 100 km level is not transparent 
for atmospheric acoustic waves — especially for long waves. Therefore, acoustic. heating is most 
effective for warming around the 100 km level rather than the 200 to 300 km levels (Figure 13), 

On the other hand, if the transmissivity of sound waves due to thermal gradient in the atmosphere 
is disregarded, most heating takes place around the 300km level and its rate exceeds 10 °C/day in 
winter. 

Although acoustic power from the polar night jet stream is more than two orders of magnitude 
larger than that of the tropospheric disturbance, which is of the order of 10' 2 ergs/cm 2 sec (Refer- 
ence 14), acoustic heating around the mesopause is not sufficient to compensate for the cooling rate 
(about 10 °C/day) of these levels in winter. 

In this respect, another mechanism of dynamical heating such as energy transfer by means of 
internal gravity waves (Reference 8) might be effective for the warming of the mesosphere observed 
by Nordberg and Smith (Reference 6). 

If the wind velocity of the jet stream exceeds 200 m/sec for several days, the layer around the 
mesopause can be warmed significantly even by acoustic waves. 
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Appendix A 
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Notation 

Acceleration of gravity, g = 98 cm/sec 2 . 

Ratio of specific heats of air, C p/ C v = 1.4, where C p and C v are the specific heat of air 
at constant pressure and that at constant volume, respectively. 

Gas constant of air, 2.87 x 10 6 ergs/ mol °C. 

Joule's mechanical equivalent of heat, 4.185 x 10 7 ergs/ cal. 

Resultant of all external forces (dyne/gm). 

Velocity vector, where u is horizontal (southward) and w is vertical (upward) component 
of air motion in cm/sec. 

The divergence of velocity per sec, i.e., 

yu yco 

v — 4- - 1 

A yx y z ' 

where x and z are horizontal (southward) and vertical (upward) coordinate, respectively. 

Small departure from static value of pressure, density and temperature, function of x 
and z. 

Static pressure, density and absolute temperature, which are the functions of ^ only in 
dyne/cm 2 , gm/cm 3 , and °K, respectively. 

Total pressure, density and absolute temperature, P = P 0 +P, p = p Q + p and T = T 0 +T. 
Heat conductivity of air (cal/cm sec°K). 

Horizontal wave length of pressure wave (km). 

Period of pressure wave (sec). 

Horizontal wave number, corresponding to k = 2tt/2 per km. 

Angular frequency of pressure wave, corresponding to r, 2n/2 per sec. 

Velocity of sound wave in the atmosphere (cm/sec). 
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H Scale height of isothermal atmosphere, H - RT 0/ g. 

D/ Dt The Euler ian derivative, y/yt + v * V. 

/x Viscosity of air, in dyne- sec/cm 2 i.e. poise (gm/cm-sec). 
v Kinematic viscosity of air, v - p/p 0 (cm 2 /sec). 
a 2 Thermal conductivity of air, a 2 = k/C p (cm 2 /sec). 
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